The Annals of Applied Probability 

2004, Vol. 14, No. 2, 612-650 

DOI: 10.1214/105051604000000044 

© Institute of Mathematical Statistics, 2004 
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A network belongs to the monotone separable class if its state 
variables are homogeneous and monotone functions of the epochs of 
the arrival process. This framework, which was first introduced to de- 
rive the stability region for stochastic networks with stationary and 
ergodic driving sequences, is revisited. It contains several classical 
queueing network models, including generalized Jackson networks, 
max-plus networks, polling systems, multiserver queues, and various 
classes of stochastic Petri nets. Our purpose is the analysis of the tails 
of the stationary state variables in the particular case of i.i.d. driv- 
ing sequences. For this, we establish general comparison relationships 
between networks of this class and the GI / GI /l/oo queue. We first 
use this to show that two classical results of the asymptotic theory 
for GI I GI /l/oo queues can be directly extended to this framework. 
The first one concerns the existence of moments for the stationary 
state variables. We establish that for all a > 1, the (a + l)-moment 
condition for service times is necessary and sufficient for the existence 
of the a-moment for the stationary maximal dater (typically the time 
to empty the network when stopping further arrivals) in any network 
of this class. The second one is a direct extension of Veraverbeke's 
tail asymptotic for the stationary waiting times in the GI / GI /l/oo 
queue. We show that under subexponential assumptions for service 
times, the stationary maximal dater in any such network has tail 
asymptotics which can be bounded from below and from above by a 
multiple of the integrated tails of service times. In general, the upper 
and the lower bounds do not coincide. Nevertheless, exact asymp- 
totics can be obtained along the same lines for various special cases 
of networks, providing direct extensions of Veraverbeke's tail asymp- 
totic for the stationary waiting times in the GI j GI j l/oo queue. We 
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exemplify this on tandem queues (maximal daters and delays in sta- 
tions) as well as on multiserver queues. 

1. Introduction. We show in the present paper that properties which 
have been known for a long time for the tail asymptotics of isolated single 
server queues can be extended to the class of stochastic networks which are 
monotone and separable. This class, which was introduced in [6], contains 
several classical queueing network models like generalized Jackson networks, 
max-plus networks, polling systems and multiserver queues. This is also 
related to the class of topical (monotone and nonexpansive) maps of [18]. 

Section 2 summarizes the definition and main results that are known on 
this class of networks, and in particular the ergodic theorems that allow one 
to determine their stability region. The notion of maximal dater is recalled. 
In a generalized Jackson network, the maximal dater is the time to empty 
the network when stopping further arrivals. In a G/G/l queue, this is just 
workload. In a FIFO tandem queue, this is end-to-end delay. 

Section 3 focuses on the proof of the moment theorem. The assumptions 
that are needed here are limited to independence. We establish the following 
generalization of the classical GI / GI /l queue moment theorem, which seems 
to be new within this setting: for all a > 1, the (a + l)-moment condition 
for service times in any monotone and separable network is necessary and 
sufficient for the existence of the a-moment for the stationary maximal dater. 

The subexponential tail asymptotic theorems are given in Sections 4 and 5. 
For surveys on the state of the art for this kind of asymptotics, see [20]. 

Section 4 gives generic upper and lower bounds which hold for all subex- 
ponential monotone separable networks and which only differ in the multi- 
plicative constants. 

Section 5 elaborates on the bounds established in Section 4. A corollary 
of Veraverbeke's theorem already proved in, for example, [1] and [2] states 
that, in the GI / GI /I queue, large workloads occur on a typical event where 
a single large service time has taken place in the distant past, and all other 
service time are close to their mean. The main new result within our setting 
is Theorem 8 which extends the notion of typical event to subexponential 
monotone separable networks; large maximal daters occur when a single 
large service time has taken place in one of the stations and all other service 
time are close to their mean. 

To the best of our knowledge, among the various classes of networks 
listed above, exact asymptotics are only known for irreducible max-plus 
networks [10]. The aim of Section 6 is to illustrate how the typical event 
theorem can be exploited to solve open questions on the exact asymptotics 
of other monotone separable networks. This is done for tandem queues in 
Section 6.1.1 and for multiserver queues in Section 6.2. 
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A first natural question is whether such asymptotics can be obtained for 
the maximal daters of all subexponential monotone separable networks. We 
have no general answer to this question yet. However, the choice of the 
tandem queue example to illustrate the potential use of the method was 
made on purpose; a tandem queue is both a generalized Jackson network 
and a reducible max-plus network. Exact asymptotics can be obtained along 
the same lines for the maximal daters of generalized Jackson networks and of 
reducible max-plus networks. These exact asymptotics require a lot of extra 
technical work, which go beyond the scope of the present paper and will be 
the object of two companion papers [7] and [8]. The exact asymptotics in 
polling systems is under investigation, too. 

A second interesting question is whether such asymptotics can be ex- 
tended to other characteristics than maximal daters. As it was shown in, for 
example, [21] maximal daters and individual waiting times may have funda- 
mentally different asymptotics. This question is addressed in Section 6.1.2 
where we show how to use the typical event theorem for monotone separable 
networks in order to derive the exact asymptotics for the stationary waiting 
or response times in individual queues of the tandem queue example. 

2. Basic results on the monotone-separable networks. 

2.1. Framework. Consider a stochastic network described by the follow- 
ing framework. 

1. The network has a single input point process N, with points for 
all m < n £ N, let iVr min i be the [m,n] restriction of N, namely the point 
process with points \T{\ m <i< n . 

2. The network has a.s. finite activity for all finite restrictions of iV; for all 
m < n E N, let X\ m ^ (N) be the time of the last activity in the network, 
when this one starts empty and is fed by iV[ m n ]. We assume that for all 
finite m and n as above, Xt m>n \ is finite. 

We assume that there exists a set of functions {/;}, fr-Mf x K l — ► K, such 
that 

(1) X [m,n]( N ) = fn-m+l{(Ti,Cl), m<l<n}, 

for all n,m and N, where the sequence {Cn} is that describing service times 
and routing decisions. 

We say that a network described as above is monotone separable if the 
functions f n are such that the following properties hold for all N. 

1 (Causality). For all m<n, 

X[ m ,n]{N) > T n . 
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2 (External monotonicity) . For all m<n, 

X[ m ,n]( N ') > X [m,n](N), 

whenever N' = f {T^} is such that T' n > T n for all n, a property which we 
will write N' > N for short. 

3 (Homogeneity) . For all c G K and for all m<n, 

X [mM (c + N)=X [mM (N) + c. 

4 (Separability). For all m < Z < n, if Xr m n(iV) < Tj+i, then 

X [mM (N)=X [l+lM (N). 

Remark 1. Single-server queues, tandem queues, and generalized Jack- 
son networks satisfy properties 1-4 above (see [5] and [6] for details). 

Remark 2 . Using the terminology of the literature on idempotency (see, 
e.g., [18]), the monotone-separable framework can be rephrased in terms 
of so-called topical forms. Indeed, for all m < n, Xt m>n -\ can be seen as 
function of the bi- infinite vector (. . . ,T_2,T_i,To,Ti,T2, . . .) of M z . Since 
Xr m>n ] : M z — > K is monotone and homogeneous, according to this terminol- 
ogy, the family X\ m)V \ , — oo < m < n < oo, is a family of topical forms on M z . 
The link between these forms is established via the separability assumption, 
which allows one to study the asymptotic forms ^(_oo,n] > which are the main 
objects of interest. Of particular interest to us are the statistical properties 
(moments, tail behavior, etc.) or the projective properties of the sequence 

(TnjX^ac^) G M 2 . 

2.2. Maximal daters. By definition, the [m, n] maximal dater is 

Z[m,n] ( N ) = ^[- m ,n] (N) ~T n = X[ m>n ] (N - T n ). 

Note that Z[ m>n \(N) is a function of {(i} m <i< n and {Ti}m<l<n-i only, where 
r n = T n+ i —T n . In particular, Z n (N) = Zy n ^ (N) is not a function of {17} 

Lemma 1 (Internal monotonicity of X and Z). Under the above con- 
ditions, the variables Xr min i and Z^ m ^ satisfy the internal monotonicity 
property; for all N, 

X [m-l,n\ (N) > X[ m ^(N), 

Z[ m -i, n ](N) > Z M (N), m < n. 

In particular, the sequence {Z^_ n ^(N)} is nondecreasing in n. Put 
^ = ^(-oo,o] = !™ z [-n,o](X) < 00. 
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Lemma 2 (Subadditive property of Z). Under the above conditions, 
{Zr m n ]} satisfies the following subadditive property: for all m < I < n, for 
all N, 

(N)<Z [m>l] (N) + Z [l+1M (N). 

2.3. Stochastic assumptions and main stability results. Assume the vari- 
ables {r n ,Cn} are random variables denned on a common probability space 
(fi, T, P,#), where 9 is an ergodic, measure-preserving shift transformation, 
such that (r n ,£ n ) o 9 = (-r n+ i,£ n+ i). The following integrability assumptions 
are also assumed to hold: 

. def ^ _ 1 dcf n 

hiT n = \ = a < oo, kjZ n < oo. 
We summarize the main results of [6]. 

Lemma 3 (0-1 law). Under the foregoing ergodic assumptions, either 
Z = oo a.s. or Z < oo a.s. 



The network is stable if Z < oo a.s. and unstable otherwise. 
Denote by Q = {T^} the degenerate input process with T' n = a.s. for 
all n. 



Lemma 4. Under the foregoing ergodic assumptions, there exists a non- 
negative constant 7(0) such that 

lim *[-„,-!] (Q) = lim vz M {Q) = ^ 

n— >oo 77, n— >oo n 

The main result on the stability region is the following theorem. 



Theorem 1. IfXj(0) < I, then Z < 00 a.s. If Z < 00 a.s., then XjiO) < 

1. 



2.4. Further assumptions. Most of the new results of the present paper 
will be obtained under the following independence assumption. 

(IA). The sequences {( n } and {r n } are mutually independent and each of 
them consists of i.i.d. random variables. 

For certain results, we shall make the following additional assumption. 

(AA). For all i, 

(2) Z i = Z [i>i] =Yf 1) + --- + Yf r \ 
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(i) 

where the r.v.'s are nonnegative, independent of interarrival 

times, and such that the sequence of random vectors (Y^ , . . . , Y^) 
is i.i.d; general dependences between the components of the vector 
(Y± , . . . , Y^) are allowed. In addition, 

o 

(3) Z [nfi] (Q)> max J2 Y P a - s - 

7=1,. ...r f — ' 

i=n 

2.5. Upper and lower bound G/G/l/oo queues. The results of this sec- 
tion are new. We assume stability, namely 7(0) < a. We pick an integer L > 1 
such that 

(4) EZ [ _ Li _ 1] (Q)<La, 

which is possible in view of Lemma 4. Without loss of generality, one can 
assume To = 0. 

To the input process N, we associate the following lower and upper bound 
processes: N~ = {T~}, where, for all k and n in Z such that n = (k — 1)L + 
1,..., 

kL, T~ = Tr^tyi, and similarly, N + = where = T^l if n = [k — 

1)L + 1, . . . , kL. Then for all n, 

(5) X { _ nfi] (N-) < X { _ nfi] (N) < X h „ )0] (AT+) = Z hns0] (N+). 

In other words, both upper and lower bound processes have batch arrivals 
(of size L). 

Note that if (IA) holds, the r.v.'s Z^ nfi ](N~) = X^ nfi ](N~) - T_ L and 
Z\-n,o]{N + ) have the same distribution and that the r.v.'s Zr n> Q\(N~) and 
are independent. 

2.5.1. Upper bound queue. The next lemma, which establishes a first 
connection between monotone-separable networks and the G/G/l/oo queue, 
directly follows from the monotonicity and the separability assumptions. 

Lemma 5. Assume Tq = 0. For any m < n < 0, 

z [ m ,o](N) < Z [nfi] (N) + max(0, Z {m ^^{N) - r n _i). 

Put Z n = Z[ n _ n ](iV). Then the sequence {Z n } does not depend on N and 
forms a stationary and ergodic sequence. 

COROLLARY 1. Assume T = 0. For any m < 0, 

/ \ 

Z[m,o] = z [m,o] (N) < max \Y. Z i~ T A 
m < k <°\ l=k i=k+1 J 

with the convention J2i = 0. 
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The main weakness of this upper bound comes from the fact that the 
corresponding queue may be unstable whereas the initial network is stable. 
This is taken care of by the upper bound described below. 

COROLLARY 2. The stationary maximal dater Z = Z/^^ (N) is bounded 

from above by the stationary response time R in the G/G/l/oo queue with 
service times 

(6) S n = ^[i( n -l)+l,in] (Q) 

and interarrival times r n = Ti n — T L ^ n _^, where L is the integer defined 
in (4). Since b = Es n < Ef„ = La, this queue is stable. 

PROOF. We have 

< supZ[_ fcL+10 ](iV + ) < sup max s + - r j+1 ) 

fe>0 A:>0 ~k<i<0 \ 



J=l 



= S + SUP V ( Si - T i+ i) = R. 

k >°i=-k 

In these relations, (5) was used to derive the first inequality, Corollary 1 was 
used in the last inequality; we also used the fact that 

Z [L(n-l)+l,Ln]{N + ) = Z[L(n- 1)+1, Ln] (Q) 

and the convention X)o 1 = 0- D 

The queue of Corollary 2 will be referred to as the L-upper-bound G/G/l/oo 
queue associated with the network. 

Note that when (IA) holds, this queue is a GI / GI /l/oo queue. In this 
case, R = W + so, where W is a stationary waiting time and W and sq are 
independent. 

Notice that under (AA), 

Ln r Ln 

(7) max r £ Y^<s n <J2 £ ^ a.s. 

3 '"'' r i=L(n-l)+l j=li=L(n-l)+l 

where the second inequality follows from the subadditive property of Z. 

The following result does not require (AA) and holds for all monotone 
separable networks such that the sequence {Zi} is i.i.d. 

We say that a nonnegative r.v. X is light tailed if there exists a positive 
number c such that Eexp(cX) is finite. 
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Corollary 3. If Zq is light tailed and A7(0) < 1, then Z^^^j * s light 
tailed too. 

Proof. From Corollary 2, it is enough to prove that the response time 
R in the L-upper-bound queue is light tailed. From well-known results on the 
G / GI /l/oo queue, in the stable case, the stationary response times R are 
light tailed when the service times are light tailed. But from the subadditive 
inequality, we have 

L 

(8) S! = Z [hL] (Q)<J2Zi, 

i=i 

which proves that s± is light tailed if Zq is. □ 

2.5.2. Lower bound fork-join queue. The following result is immediate. 
Lemma 6. Under Condition (AA), 

/ o -i \ 

(9) Z^oj > R = max sup £ - ]T < 

J— l,—,r n <o \ 



The queue with service times {Y^ } and interarrival times {r^} will be 
referred to as the j'-lower-bound G/G/l/oo queue associated with the net- 
work. Let By' denote the stationary response time in this queue: 

i^ = su P (£yW-£4 

™<0 \ n n J 

Then the lower bound R defined in (9) is the stationary response time in 

the r-dimensional fork-join queue with service times {Y^ '}, j = 1, . . . , r and 
interarrival times {t^}. 

2.6. Examples. 

2.6.1. Tandem queues. Consider a stable G/G/l/oo — > -/G/l/oo tan- 
dem queue. Denote by {o - !^} the sequence of service times in station i = 1, 2 
and {r n } the sequence of interarrival times at the first station. Put W' = 
Eo-«, a = Er and p^=b^/a < 1. We have 7(0) = max^ 1 ), b^). 

Tandem queues fall in the class of open Jackson networks, and in the 
class of open max-plus systems which both belong to the class of monotone 
separable networks (see below) . We have the following representation for the 
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maximal dater (see, e.g., [10]) 

(10) Z { _ nfi] = sup sup (t^+E^-^o-Tp)], 

-n<p<0p<q<0 \ m=p m=q J 

(11) Z = = sup sup ( °m + E CT - 2) " ( T o " Tp)) . 

p<0p<q<0\ m=p m=q J 

Assumption (IA) is satisfied if the sequences {r n } and {£„ = (an\ cr^)} 
are i.i.d. and mutually independent (we may allow a dependence between 
an"* and <Jn^ ) ■ Assumption ( AA) is also satisfied here with r = 2 and Yn = 
a^, i = l,2. 

The maximal dater with index n is the sojourn time of customer n in the 
network, namely the time which elapses between its arrival in station 1 and 
its departure from station 2. 

As for the L-upper-bound queue associated with this network, the expres- 
sion for s n is here 

( 12 ) *n = ™f< X L ( E vfn-DL+i + E ^IjL+il " 

— \i=l i=j / 



2.6.2. Multiserver queues. Let 

Ty n = (^w,...,^ m )) 

be the Kiefer-Wolfowitz workload vector in the GI / GI /m/oo queue with 
inter arrival times {r n } and service times {<r n }. Here n is the customer index 

(i) 

and Wk>, i = l,... ,m, are the workloads of the servers at the nth arrival 
time, arranged in nondecreasing order. More precisely, we assume Wq = 
(0,...,0) and 

(13) W n+1 = Yi{W n + ei a n - ir n )+ 

for i > 0, where ei = (1, 0, . . . , 0) and i = (1, 1, . . . , 1) are m-dimensional vec- 
tors and the operator R permutes the components of a vector in nonde- 
creasing order. For a multiserver queue, 7(0) = ~Eao/m. 

Assumption (IA) is satisfied under the assumption that the service times 
are i.i.d. Assumption (AA) is not satisfied here. 

The maximal dater associated with customer n is the time which elapses 
between its arrival and the time when all customers still present at its arrival 
time have left the system (including customer n), 



Z [0in] =max(W« +a n ,wi m) )- 
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2.6.3. Generalized Jackson networks. Consider a generalized Jackson net- 
work with r stations. We denote by: 

1- {<7n^} the i.i.d. sequence of service times in station k. 

2. {fi^} the i.i.d. sequence of routing decisions from station i; with values 
in the set {1, . . . , r}. 

3. {fin} the i.i.d. sequence of routing decisions for the input process; also 

with values in the set {1, . . . , r, r + 1}, where fin = r + 1 means that a 
customer taking nth service at station % leaves then the network. 

4. {r n } the i.i.d. sequence of interarrival time. 

Under these assumptions, both (I A) and (A A) are satisfied. We have 

(i4) n w = E^ 

l 

with v{j) the total number of visits of customer 1 (the customer arriving at 
time Ti) to station j in the [1, 1] restriction of the network, namely when this 
customer is the only one to enter the network. The random variables v(J), 
j = 1, . . . , r, are obtained from the sequences of routing decisions (see [5]). 

In this case Zi n> Q-\ is the time which elapses between the arrival of cus- 
tomer and the time when all customers have left the system, given that 
arrivals are stopped after To. 

2.6.4. Max-plus networks. The class of open max-plus networks also falls 
in this framework (see, e.g., [6]). A typical example of this class is that of 
tandem queues. Tandem queues form a reducible open max-plus network. 
For examples of irreducible networks of this class, see [10]. 

3. Integrability of stationary maximal daters. We assume (IA) and sta- 
bility, namely A7(0) < 1. 

Let W denote the stationary waiting time in the L-upper-bound GI/GI/l/oo 
queue of the network. The following result is well known. 

Lemma 7. For any a > 1, EM^ a_1 is finite if and only i/EsQ is finite. 

Therefore, R = W + s is such that Ei? Q_1 is finite if and only if Esq is 
finite. 

Corollary 4. IfEZfi < oo, then EZ"'^ , < oo. 
PROOF. We have 

o 

so < Zi - 

-L+l 
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Therefore if EZft < oo, then Esq is finite. Thus, E^"" 1 and E^ " 1 are 
finite, too. We conclude the proof by using the bound Z/^q] < R (see the 
proof of Corollary 2). □ 

Under condition (AA), EZq is finite if and only if for all j, E[(Zo(j)) a ] is 
finite. The following theorem is then an immediate consequence of Lemmas 
6 and 7. 

Theorem 2. Under assumptions (IA) and (AA), j/EfZ^o ]] * s finite, 
so is EZ$. 

Examples. All results are given under the assumption that the system 
under consideration is stable. 

1. Tandem queues. The system response time has a moment of order a — 1 
iff the service times in both stations admit a moment of order a. 

2. Multiserver queues. In steady state, the time to empty the system has a 
moment of order a — 1 if the service times admit a moment of order a. 

3. Generalized Jackson networks. The stationary maximal dater has a mo- 
ment of order a — 1 iff all service times have moments of order a. Since 
the stationary maximal dater is not less than the residual workload at 
any station, we get if all service times have moments of order a, then 
the stationary residual workloads and the stationary queue lengths at 
all stations have moments of order a — 1. Since the number of customer 
services has a geometrical tail, one also deduces from this that stationary 
sojourn times also have moments of order a — 1. 

4. Bounds for sub exponential tail asymptotics. 

4.1. Assumptions and notation. Here and later in the paper, for strictly 
positive functions / and g, the equivalence f{x) ~ dg{x) with d > means 
f{x)/g{x) — > d as x — > oo. This equivalence may also be rewritten as f(x) = 
dg(x)(l + o(l)) = dg(x) + o(g(x)) = dg(x) + o(f(x)), where o(l) is a function 
which tends to as x tends to oo, and o{g{x)) is a function such that 
o{g{x))/g{x) — > as x — > oo. By convention, the equivalence f(x) ~ dg(x) 
with d = means f(x) = o{g{x)). We will also use the following notation: 

1. f{x) = @(g(x)) to mean lira sup f (x) / g(x) < oo and liminf f(x)/g(x) > 0, 

2. f(x) = 0{g{xj) to mean limsup / '(x) / 'g(x) < oo. 

4.1.1. Tails. Let £ be a nonnegative r.v. with distribution function F 
such that P(£ > x) = 1 — F(x) = F(x) > for all x. Let £1,^2 be independent 
copies of £. 
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Definition 1. £ has a heavy-tailed distribution (HT), if, for any c> 0, 

("OO 

Eexp(c£) = / exp(cx) dF(x) = oo. 
Jo 

Definition 2. £ has a long-tailed distribution (LT), if, for any y > 0, 
F(x + y) ~ F(x) as oo. 

Any LT distribution is HT. 

Definition 3. £ has a subexponential distribution (SE), if 
P(&+6>z)~2F(a;) asx^oo. 

Any SE distribution is LT. For basic properties of subexponential distri- 
butions, see, for example, [13]. 

4.1.2. Network assumptions. Consider a distribution function F on K + 
such that the following hold: 

(a) F is subexponential, with finite first moment M = J °° F(u) du, where 
F(u) = 1 — F(u) denotes the tail of F. 

(b) The integrated tail distribution F s , 

F s (x) = l-mmll,J™F(u)du\ = 1-F s (x), 
is subexponential. 

Here are a few properties satisfied by F that will be needed later on and 
that follow from the fact that F s is long tailed. 
When x — > oo, 

(15) F{x) = o(F s (x)). 

As a corollary, there exists a nondecreasing integer-valued function N x — > oo 
and such that, for all finite real numbers b, 

N x _ 

(16) ^F(x + n6) = o(F"(x)), X — > oo. 

n=0 

In particular, 

(17) N x F(x) = o(F s (x)). 

Such a c.d.f. F being given, consider a monotone separable network satis- 
fying (IA) and (AA) and such that the following equivalence holds when x 
tends to oo: 
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(c) For all j = 1, ... , r, 

p(y 1 (i) > x) ~ ^T(x) 

with X)j = d > 0. 

For a monotone separable network, the three assumptions (a)-(c) will be 
referred to as (SE). Under (SE), the following holds: 

(18) / P(y i (j) > y) dy ~ d^F 8 (x) as x — ► oo. 

Assumption (H). We also introduce the following assumption (H): 

p ( £ y i (j) > x l ~ p ( y i (j) > x 

1 1 

Note that the very last equivalence follows from (SE). Assumption (H) is, 
for instance, satisfied in the particular case when the random variables 
are mutually independent; in Section A. 2, we give sufficient conditions for 
(H) to hold that go beyond this particular case. 

Take any 1 < i±, %i < r, i\ ^ %<i. Since 



x 



P^maxy^ > xj <Y, P ( Y i j) >x)- P(y x (il) > x,Y} l2) > 
we deduce from (19) that 

(20) P(Y x (il) >x,Yl i2) >x) = o(F{x)). 

Remark 3. In what follows, we will not need i.i.d. assumptions on the 
interarrival times {r n }. As it follows from Theorem 14, the results we prove 
will hold also in the more general situation when these variables satisfy the 
following three conditions: 

1. {r n } forms a stationary ergodic sequence with a finite positive mean 
Eti = a. 

2. {r n } is independent of {Yn\j = 1, . . . ,r}. 

3. For all a < a, 

P sup na — Tj > x = o( F s (x)). 
\n>o\ i= _ n J J 

(See [4] for the proof in the single-server queue case.) 
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4.2. Tail asymptotics for the supremum of a random walk with subexpo- 
nential increments. We now remind the well-known result from [14] and [22] . 
We use negative indices in order to link the result with queueing applica- 
tions. 

Theorem 3. Let {£ n } be an i.i.d. sequence with negative mean E£i = 
—a, Sq = 0, S- n = J2iC-i an d S = sup n>0 S'_ n . Assume that there exists a 
distribution function F on [0, oo) such that F s is subexponential and P(£i > 
x) ~ dF{x) with d> as x — > oo. Then, as x — > oo, 

P(S>x) = (l + o(l))-F s {x). 

a 

In particular, consider a GI / GI/l/oo queue with i.i.d. service times {a n } 
(with mean b) and i.i.d. interarrival times {r n } (with mean a> b) and put 
£n = o~ n — T n . Assume that P(o"i > x) ~ dF(x), with F as above. Then the 
stationary waiting time W and the stationary response time R of customer 
are such that 

P(R > x) ~ P(W > x) = (1 + o(l)) -^—F s (x). 

a — b 

In particular, if the distribution function of a is F , then P(i? > x) ~ P(W > 
x) = (l + o(l)) 5 i 5 F(x). 

The following lower bound is also known (see, e.g., [4]) and was obtained 
by the use of the strong law of large numbers (SLLN): 

Theorem 4. Consider a G/G/l/oo queue with i.i.d. service times {cr n } 
(with mean b) and independent stationary ergodic interarrival times {r n } 
(with mean a> b). Assume that P(ui > x) ~ dF{x), where d > and where 
the integrated distribution F s is long tailed. Then 

P(R >x)> P(W >x)> -^—F s (x) + o(F s (x)). 

a — b 

4.3. Bounds. 

4.3.1. Upper bound. Let Z = 2y_oo o] and let L be the integer defined in 
Section 2.5 and let s be the service time in the associated L-upper-bound 
GI I GI j\ I oo queue. 

Put 6 = Es and note that Ef = La. Then p = ^ = \j(0)(l + 6) < 1 where 
5 may be chosen as small as possible. We deduce from (7) and (H) that 

P(si >x)~dLF(x). 
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Thus, from Theorem 3, 



l r°° 

P(R>x)~- ? / P(s>y)dy 

La — b Jx 



La — b Jx 

1 r<*> , d 



f°° d — 

/ dLF(y)dy = ^— F s (x) 

Jx a — b/L 



La — bJx a — b/L 

Here b/L — ► 7(0) as L — > 00. We have proved the following. 

Theorem 5. Under the (IA), (AA), (SE) and (H) assumptions, 
(21) limsup — — s - < lim lim — L=g 



F (x) i? A ( x ) a -7(0)' 

Remark 4. The assumptions of Theorem 5 bear on the random vari- 
es 

ables Y{". These can be weakened by considering conditions on the random 
variables Z n = Zi n>n i as follows: If the random variables Z n are i.i.d. with dis- 
tribution G such that both G and G s are subexponential, and if the random 
variables r n are i.i.d. and independent of the {Z n } sequence, then 

, , , P(Z>x) 1 

(22) limsup — s < 



x^oo G (x) a - 7(0) 
The proof of this is based on Corollary 2 and on coupling arguments. 

4.3.2. Lower bound. From (9), 

Z = Z(_ OO)0] > R = max sup ( J] - ^ n \ = max . 
Then from Theorem 3, 



3 n<0\ i=n i=n I J 



P(R®>x)~ ( -, F S (x), 
v ' 0-6W 

with &&') = EF/ j) . Note that 

^P(i? (j) >x) >P^maxi? (j) > 

(23) ^ 

> J2 p ( rU) > x ) - Yl p ( R{il) > R(i2) > x )- 

Since, for any i±^i2, 

(24) P(R^ >x,R {i2) >x) = o(F s (x)) 
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(see Section A.l for the proof), we get 

P(R>x) =Y,P{R U) >x) +o(F s (x)). 

j 

Thus, the following theorem holds. 

Theorem 6. Under Assumptions (IA), (AA), (SE) and (H), 

(25) Kmi ,f > lim ? = ± -#L. 

y ' *-+°o F (x) - «-»oo F (x) 4^ x a-tf3) 

Remark 5. The asymptotics for the lower and upper bounds are the 
same up to multiplicative constants. So Theorems 5 and 6 imply P(Z > x) = 
G(F s (x)). 

In the single-server isolated queue case, 7(0) = b = M, b = Lb and d = 1. 
Therefore, in this case the upper and lower bounds coincide. 

4.4. Examples. 

4.4.1. Tandem queues. The definitions and notation are those of Sec- 
tion 2.6.1. We assume that 

(26) Fi(x) = P(<7« > x) ~ d®F(x), 

that d = d^ > and that both F and F s are sub exponential. Assump- 

tion (H) is valid if we assume in addition that o~n^ and <ri are independent. 

Denote by Z the stationary sojourn time in the network. We look for the 
asymptotic behavior of the function P(Z > x) as x — > 00. 

The lower bound (25) is 

r P(Z>x) dP> dW 
hmmf -4=i= > -7TT + 



F\x) ~ a-bi 1 ) a -6( 2 )' 
Since 7(0) =b = max(6^ 1 \ b^), the upper bound (21) reads 

P(Z>x) d 
hmsup =j— < -. 

x^oo F (x) a — b 

This upper bound was proved in [10]. 

4.4.2. Generalized Jackson networks. The definitions are notation and 
those of Section 2.6.3 and of Section A. 2 in the Appendix. We assume that 
P(crW > x) ~ l^F(x), that ~£j > and that both F and F s are subex- 
ponential. Put ir® = Eu^) and b^ = EY® = tt^Ect^. Without loss of 
generality we may assume b^ to be positive for all j. The network is stable if 
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7(0) < a, where 7(0) = b = max(6^ 1 ), . . . , b^ r ') (see, e.g., [5]). Assumption (H) 
is valid (see the example at the end of Section A. 2). 

Denote by Z the stationary maximal dater. Then, from (25), the lower 
bound for P(Z > x) is 

,. . V(Z>x) 'l®*® 
hm mi =g- - — > > -7TT- 

x ^oo F *( x j -Z-J a _ h (j) 

and the upper bound (21) reads 

P(Z>x) £Z^V^ 
hmsup a < — . 

x->oo F (x) a — b 

4.4.3. Max-plus networks. Similar bounds were studied within the frame- 
work of open, irreducible max-plus networks in [10] . As in the single-server 
isolated queue case (which is an instance of such networks), the upper and 
lower bounds coincide, which yields the exact asymptotics. However, the 
exact asymptotics are not known for reducible max-plus networks, even for 
the particular case of tandem queues. 

5. Typical event of a subexponential monotone-separable network. 

5.1. Typical event of a subexponential GI / GI /l queue. This section con- 
tains qualitative indications on how rare events occur in subexponential 
GI I GI 1 1 queues in terms of asymptotic equivalences involving the so-called 
typical event. 

Results of the same nature were first stated by Anantharam in [1] in 
the regularly varying case (see Theorem 3.1 therein) and by Asmussen and 
Kliippelberg in [2]. However, we could not find the equivalences (Corol- 
lary 5) precisely needed for the extension to monotone separable networks 
(Theorem 8) in any earlier paper. Some notation and ideas of the proof of 
Corollary 5 will be used in Theorem 8. 

Consider a GI j G//I/00 queue with mean interarrival times a = Er n and 
mean service times b = E<x„, where a> b. Denote by F the distribution func- 
tion of a and assume that F satisfies (b) of (SE) (i.e., F s is subexponential), 
and let N x be the associated function defined in Section 4.1.2. Let 

n n n 

£,n = 0- n — T n , = yV-j, S'^ = ^0'_i, S n = / J £-j =S° — 5^. 

1 1 1 

Corollary 5. Let W (resp. R) denote the stationary waiting (resp. 
response) time of customer in the FIFO GI/GI/l/oo queue. For any x, 
let {K niX } be a sequence of events such that: 
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(i) For any n, the event K n ^ x and the random variable c_ n are indepen- 
dent; 

(ii) ~P(K nx ) — > 1 uniformly in n> N x as x — > oo. 
For any sequence r\ n — > 0, let 

(27) A n>x = K ntX r\{a- n >x + n(a- b + %)} and A x = (J ^4 ni2; . 

n>JV» 

Then, as x — > oo, 

P(W > x) ~ P(W > x, A,) ~ P(A,) 



(28) 



oo oo 
n=N x n=N x 



J2 ?{W>X,A n>x )~ J2 P ( A n, 



and 

(29) P(R>x)~P(W>x). 

Proof. Simple calculations using the fact that F s is long tailed show 
that, as x — > oo, 

P(A n ,x)= P(Kn,x)V(o-- n >x + n(a-b + r) n )) 

n>N x n>N x 

V F(x + n(a-b + r, n ))~-^—F s (x). 

* — * n — h 



a — b' 

n>N x 



Thus, if the sequences {K UiX } and {rj n } are such that, for all sufficiently 
large x: 

(a) the events A n ^ x are disjoint for all n > N x ; 

(b) A n>x <^{W>x} for all n>N x ; 

then 

P(T^ > x) > P(W > x, A x ) = P(A X ) 

= J2 V(W>x,A niX )= J2 P(A n , x )~-L-F s (x). 

n>N x n>N x 

Combining this with Theorem 3, we get the equivalences (28). 

We now construct two specific sequences {K n>x } and {r) n } satisfying (a), 
(b) and the assumptions of the corollary. Due to the SLLN, there exists a 
nonincreasing sequence e n — > such that ne n — ► oo and, as n — > oo, 
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Put r] n = 3e n and 



<e k VN x <k<n}n 



-r--b 



<e k+1 VN x <k<n 



k 

Clearly, the conditions of the corollary are satisfied. Since ne n > b for all 
sufficiently large n, on the event A n>x , 

W > S n > x + nr] n — (2n — l)s n — b>x. 

In addition, the events A UtX , n > N x , are disjoint if £m x < (a — b)/2. In- 
deed, on the event A HiX , we then have S n > x and S^-i = maxo<j< n -i Sj < 
maxo<j< n _i j(b — a + 2sn x ) < 0; and the events {S*_ ± < 0} n {S n > x} are 
obviously disjoint. 

Take now any other sequences {K UtX } and {fj n } satisfying the conditions 
of the corollary and denote the corresponding events by {^4 n ,x} and A x . 
Then 

\P{A X )-Y{A X )\ 

<P (J {a^ n > x + n{a-b + mm(rj n ,rj n ))}\ 

V n>N x ) 

-P |J nf v n{cr_ n >x + n(a-6 + max(?7 n ,?7 n ))} 

V n>N x J 

< F ( a -n >x + n(a- & + min(77 n , %))) 

n.>JVa- 

~ 51 P (- ftr n,xni^„, a; n{o-„ n >x + n(a-6 + max(7/ n ,r? n ))}) 

n>A^ 

< Aj. F(x + n(a-b + max(r] n ,rj n ))) 

n>N x 

n(a - b + min(r] n ,r] n ))) 



n>N x 

+ 



n>N x 

- F(x + n(a - 6 + max(?7 n , r/ n ))) 

where A x = sup n >jv.,.(P(-K^ x ) + P(-K"^ ia ,)) — ► as x — ► oo . Thus, both terms 
in the last expression are o(F s (x)), and the equivalences (28) hold for se- 
quences {K n;X } and {rj n }. 

Finally, the equivalence (29) follows from the relation R = W + ao, the in- 
dependence and the fact that the tail of W is heavier than that of a [see (15)] . 

□ 

The event A x (which will be referred to as the typical event of the subex- 
ponential GI/GI/1 queue in what follows) occurs if there is only one big 
service time and all other service times or interarrival times follow the SLLN. 
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5.2. Key equivalences for the maximal dater of a subexponential monotone- 
separable network. In this section, we consider a monotone-separable net- 
work satisfying (IA), (AA), (H) and (SE). The function N x is that associated 
with the reference distribution function F of the (SE) assumptions. 

Theorem 7. Let Z be the stationary maximal dater of some monotone 
separable network. Denoting A x the typical event of the L-upper-bound queue 
(more generally we will add a hat to indicate that a variable pertains to the 
upper bound queue), we have 

oo 

(31) P(Z > x) ~ P(Z > x, A x ) ~ Y, P ( Z > X > An,x) 

n=N x 

and 

(32) P(Z>x,A x ) = Q(F\x)). 
Also, for any random variable Z such that Z < Z a.s., 

(33) P(Z>x) = J2 P(Z>x,A n , x ) + o(F s (x)). 

n>N x 

PROOF. Since Z <R a.s., 

P(Z > x) = P(Z > x, A x ) + P(Z > x; R > x, {A x ) c ) 

= P(Z>x,A x ) + o(F s (x))= Y, P(Z>x,A rhX ) + o(F s (x)) 

n>N x 

from Corollary 5. 

From Theorems 5 and 6, P(Z > x) = Q(F (x)). Thus, 

P(Z > x,A x ) = P(Z >x)- o(F s {x)) = Q(F\x)) - o(F s (x)) = @{F~ s (x)) 

and both (31) and (32) follow. □ 

The main result of the paper concerning subexponential monotone sep- 
arable networks is the following theorem, which can be seen as a network 
extension of Corollary 5 and which gives the shape of the typical event 
creating a large maximal dater in such a network. 

Theorem 8. The assumptions are the same as in Theorem 7. Put b = 

(i) 

7(0) . For any x and for j = 1, . . . , r, let {Kn,x} be a sequence of events such 
that: 

(i) For any n and j, the event Kn}x and the random variable are 
independent; 
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(ii) For any j, P(Kn)c) — > 1 uniformly in n> N x as x — > oo. 
(7) 

For all sequences r]n , j = 1, • • • ,r, tending to 0, put 



(34) 



4& = K & n { y n > x + n(a - 6 + t#>)}, 

00 r 
4° = U 4??x and A X =\JA^. 

n=N x j=l 



Then, as x — > 00, 

r 

P(Z > s) ~ P(Z > x, A x ) ~ 51 P ( Z > ^) 



(35) 



1 

r 00 

]T £ J>(Z>z,A®). 

j=ln=N x 



Similarly, for any random variable s.t. Z < Z , 

P{Z >x) = P(Z > x, A x ) + o(F s (x)) 

(36) =j2^(Z>x,A^) + o(F s (x)) 



1 

r 00 



= £ £ P(Z>x,A^ x ) + o(F s (x)). 

j=t n=N x 

IfP(Z>x)=Q(F(x)), one can replace the last equalities by equivalences 
and delete the o(F (x)) terms in the last relation. 

The equivalences (31) and (35) will be the key relationships for the exact 
asymptotics of the examples of Section 6. They show that for the monotone 
separable network also, whenever the maximal dater is large, at most one of 
the service times is large whereas all other ones are moderate. 

Proof. We will only prove the equivalence 

r 00 

(37) ?{Z>x)~Y. E *(Z>x,4&). 

j=l n=N x 

The other equivalences in (35) may be obtained similarly. 

We start the proof with the following three reductions. 

First, it is sufficient to prove the equivalence (37) when replacing each of 
the events Kn\ by the whole probability space O. Indeed, put An} x = {X-n > 
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x + n(a — b + rjn^)}. We know from Theorem 7 that P(Z > x) = @(F s (x)). 
Suppose that 

p(z>x)~j2 E p(^>z,4&)- 

j=ln>N x 

Then 

£ e p(z>x,Ag)j 

=EE p(^>^,iS)-E E p(z>z,iS\4^). 

j=ln>N x j=ln>N x 

The result then follows from the fact that the last subtracted sum is non- 
negative and is not bigger than 

^^P(ig) )P((^)) C ) < A(x)e(F s (x)) = o(F s (x)) 

since A(x) = maxi<j< r sup n>Nx P{{Kn) x ) c ) — ► as x —> oo. 

Second, it is sufficient to consider the case r/n = for all n and j. This 
follows from the following bound where S x = maxj sup n >jy (jjn ) + : 

r 

E V(Z>x,Y^>x + n(a-b)) 

j=ln>N x 

- P(Z > x, > x + n{a - b + (??^) + )) 

r 

<E E P(^-n ) e(^ + ^(a-&)^ + ^(a-& + ^))) 

= (1 + o(l)) dF s (x) 1— ) = o(F s (*)) 

\a — b a — b + o x J 

(i) 

and a symmetrical bound for the negative part of ijn ■ Thus it is enough to 
prove the equivalence 

r 

(38) P(Z>x)~J2 E P(^>a:,l r i2>a; + n(o-6)). 

j=l n >N x 

Third, if there exists a sequence sl G (0,a — 6), sl — ► such that, for 
any L, the following equivalence takes place (where bi = b + e^): 

r 

(39) P(Z>x)~^ ]T P(Z>x,yi J n } >x + n(a-6 L )), 

j=ln>AT :I 
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then (38) holds. Indeed, take el < (a — b)/2. Then 

r 

£ J2 P{Z>x,Y®>x + n(a-b L )) 

3=1 n>N x 

- P(Z > x, yjj >x + n(a- b)) 



<Y. P{Y U }e[x + n(a-b L ),x + n(a-b)]) 

j=ln>N x 

1 1 



{l + o{l))dF"(x) 



a — bL a — b 



<^ + o{l))7^F\x). 



2e L d ■ 

(a -by 

Letting L — > oo, we derive (38) from (39). 

Before proving (39), we recall that, from conditions (SE) and (H), 



P( ?x > x) ~ P ( U (J {Y^ > £ 2 P(Y^ > x). 

\j=ll=l ) 3=\l=\ 

Since F is long tailed, we can replace the latter equivalences by 
P(s l >x)^p((j \J{ Yl U) >x + l(a-b L )}\ 

~^^P(y/ j) >x + /(a-6 L )). 

3=1 1=1 

More precisely, when denoting the event in the left-hand side by C x and the 
event in the center by D x , we get C x C D x and 

(4i) sup P(A / \c,) =o(1)i 

S/>z ^(y) 

when x — > oo, while P(C X ) = 0( -F(x)). 

We now prove (39). For any L, put N x = N x /L (more precisely the integer 
part of this ratio) and note that N x also satisfies condition (16). Take L 
sufficiently large and set b^ = Eso/L = b + sl- For the L-upper queue, one 
can take the typical event of the form 

A x= U {s- n >x + nL(a-b L )}= (J A n ^ x . 

n>N x n>N x 

From Theorem 7, 

P(Z>x)~ P(Z > x,s„ m > x + mL(a - 6 L )). 

m>N x 
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From (40) and (41), 

P(Z > x,s- m > x + mL(a — 6/,)) 

= (l + o(l)) 

x E J2Y, P ( Z>x > Y -iL + i>x + ™L(a-b L ) + Ka-b L )), 

m>N x i= ll=l 

where the uniformity in m required to obtain the term o(l) follows from the 
uniformity in y in (41). So 

r 

P(Z>x)~Y,Y, P ( Z > Y - ] n+i >x + n(a- b L )) 

j=l n >N x 
r 

~E E P{Z>x,Y®>x + n(a-b L )). 

j=l n>N x ^ 

6. Two examples of exact tail asymptotics. This section gives two il- 
lustrations of the use of Theorem 8 in order to derive exact asymptotics. 
Without loss of generality we can assume interarrival times to be constants 
and equal to a (see Section A. 3). 

6.1. Tandem queues. For tandem queue, the assumptions on the tails 
of the service times are those of Section 4.4.1. Here Yn = On , so that 
(H) trivially holds, since the service times are independent. The results are 
stated for the two-station case, but the extension to tandems (or treelike 
networks) of any dimension is immediate. 

Choose a sequence N x satisfying (16). Let Wn be the stationary wait- 

(2) 

ing time of customer n in queue j = 1,2, and rA be the interarrival time 
between the nth and (n + l)st customers to the second queue. Let Cn ] = 

<7n^ —T n = an** - a, ffi = - r,i and £ n = affl - a^li- We will also use 
the following notation: 

n n n 

§n=E^ s^=y:&, s^=y:°% 3=i,2- 

i=l i=l i=0 

The following relations hold: 

(42) = max (0, + $>), j = 1, 2, 
and 

(43) = - min (0, W« + £ W ) + a % > a£l 
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so that tP = aJjJk if WP + } > 0. In addition, 

(44) Z = w! ) 1) +a^ + wP+a^. 
Also, from (11), 

(45) Z= sup (j^a^ + S^+pa). 

-oo<p<g<0 \ m=p / 

Similarly, for any n, the stationary response time Z^^ _ n j of customer (— n) 
satisfies the relations 

%00,-n] = ^-n + + ^ + (7® 

(46) / q -n \ 

sup E ffW + ^ a$ + (p + n)a . 



-oo<p<<7<— n 



, m=p m=q 



6.1.1. End-to-end delay. In this section, we prove the following exact 
asymptotic, which refines the bounds of Section 4.4.1 (these bounds do not 
coincide in general). 

Theorem 9. Under the assumptions of Section 4.4 ; 

/ j(1) j(2) \_ 

(47) P(Z>I) .(_ + _) F ' W , 

w/iere 6 = max(&W , ) = 7(0) . 

Remark 6. As a corollary of Theorem 9 and of results from [3] and 
[16], one can easily derive sharp asymptotics for the stationary queue length 
Q = Qi + Q2 in the tandem queue. Also, the result may be easily extended 
to queues in tandem of any finite length. 

Proof of Theorem 9. Prom Theorem 8, we get 

2 00 2 

p(z>*)~£p(z>z,4^v E E p (^>^S)- 

j=l n=N x 3=1 

We have to find appropriate sequences {Knl} and {rjn^}. 
Start with j = 1. For any {Kn)x} and T]n — ► 0, 

E P(^>*,4i)< E P(^>x + n(a-6) + ?? ( 1 ))~-^-r(x). 

n> n>7V„ 



26 



F. BACCELLI AND S. FOSS 



For the lower bound, consider the events 

#S = {flfc? > n(bW - r£)),S™ > n{b^ - ,£>)} 

and choose a sequence r]n^ — ► such that P(lfn*x) — > 1 uniformly in n > N x 
as x — > oo. Then, from (45), 

P(Z > x, A« ) > P(a« + max (S**?, S<T' 2 >) - na > x, A«) 

> P(a a l + n(max (&« ,6^) - - a) > x)p(ifg) 

= (1 + o(l))P(a^ > x + n(a - 6 + r?«)), 

and the lower bound for P(Z > x,A|c ) is asymptotically equivalent to the 
upper one. 

Consider j = 2. The lower bound 

P(Z > x,4 2) ) > P(^o (2) > 4 2) ) = (1 + o(l))-^F s (x) 

follows from Theorem 4. 
For the upper bound, put 

U n = sup sup ( °in ] + Si a ' 2) +Pa) 

-oo<p<0max(p,-n)<(j<0 \ m =p J 

and note that, from (45) and (46), 

Z < max (Z(_ ao _„_!] + 5^ CT ' 2) - na, £7 n ) 

= max (Z^oq _ n _i] + + a^l - na, U n ), 

where the random vector (Zr^ ,-n-i]> U n , ) is independent of cr_„,. 

Since C/ n < Z a.s., P(U n < x) — > 1 uniformly in n as i -> oo. Since the 
distribution of Z(_ oc „ n _ 1 ] does not depend on n, Z^^^^i/n — > in 

probability. Due to the SLLN, S^J 2 ^ /n — > b^ a.s. Therefore, there exists 
a sequence e n J. 0, ne n — ► oo such that 

P(£7 n < x, Z ( _ 0O| _ n _ 1] < ne n , S<£? < n(6^ + e n )) -► 1 

(2) 

uniformly in n > iV x as x — > oo. Denote the latter event by Kn,x and recall 

(2) (2) 

that it is independent of cr_„. Put ijn = — 2e n . Then 

P(Z>x,4 2 ).)<P(a^ + n(6( 2 ) — a) + 2ne n > x, A^,) 

= P(a^ n >x + n(a-b^ + ^),Ki% 

= (1 + o(l))P(^ > x + n(a - + r^)), 
and the desired asymptotics follow. □ 
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6.1.2. Delay at the second queue. In this section, we focus on the asymp- 
totics for the stationary waiting time = of customer at the 

second queue. The assumptions are the same as in Section 6.1. 

Results on the matter were obtained by Huang and Sigman in [19] in the 
case where the tail of heavier than that of a^ l \ The results of the 

present section are more general in that such an assumption is not required. 

First, let us see how the results of [19] follow from what we have here. 
Under the assumptions of Section 4.4.1, we get from (36) of Theorem 8 (for 
Z = < Z) that 

(48) P(W {2) >x) = P(W {2) >x,A^) + P(W {2) >x,A^) + o{F s {x)). 
Then 

,7(2) 

(49) P(W® > x , A f) = — W)F s (x) + o{F s {x)). 

The lower bound follows from Theorem 4 and the upper one from the in- 
equality P(W^ > x, A {2) ) < P(Z > x, A (2) ) and from the part j = 2 of the 
proof of Theorem 9. For the first term in the right-hand side of (48), we have 

r/(l) _ 

V<P{W {2) >x,A { i ) )<P{Z>x,A i £ ) ) = —-F s {x)+o{F s {x)) 

a — b 

from the proof of Theorem 9. Thus, if = < S 2 \ then 

,7(2) 

(50) P(W®>x) tm fS ( x ) 

which is the result of [19]. 

We now successively consider the three cases 

foW > 6(2), b W = b W ail d 

6(D < 6 (2). 

Case 6( 1 ) > 6( 2 ) . For the following theorem, we do not need any assump- 
tion on the tail of cr( 1 ). In fact, we do not even need to assume that F is 
subexponential, the only required assumption being that F s is subexponen- 
tial. 

Theorem 10. Assume a > > b^ and P(<j( 2 ) > x) ~ d^F(x) where 
c?( 2 ) > 0, and the integrated tail distribution F s is subexponential. Then, as 
x — > CO, 

(51) P{W W>x)~-?-Lf s (x). 
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Proof. We already established the right lower bound in (49). Thus, it 
is enough to derive an upper bound which is asymptotically equivalent to 
the lower one. 

For this, we use the notation from the beginning of Section 6.1. Since 

(2) (2) 

a I and t( are independent and F s is long tailed, as x — > oo, 
(52) / P(ef } > t) dt ~ / P( ( r i (2) > t) dt ~ d^F 5 ^). 

Put 5( 2 ) = sup n > Si 2) and 5 = sup n > S n . Since r/ 2) > a\% = < 
S a.s. Since > b^ 2 \ S < oo a.s. So, we have 

P(S( 2 ) >x)= P(S (2) > x, S > x) 

(53) 

<]TP(S (2) >x,^ n >x + nc) + o(F s (x)), 

n 

where (53) follows from Corollary 5 which implies that 

{S>x}= U-te-n >x + nc}U B x where ¥{B X ) = o{F s {x)). 

n 

Set c = a - b^ and c = - 6 (2) . For all e G (0,2), i? > and n, define 
the event 

A.,*,* = {sf } < R - i(c - e), Si < R - i(c - e), i = 1, 2, . . . ,n - 1; 

5 n+J - - S n < R- j(c-e), i = 1, 2, . . . }. 

By the SLLN, for any e > 0, there exists R > such that, for any n = 1, 2, . . . , 
P(-Cn,e,i?) > 1 — e. From (53), we have 

P(5< 2 > >x)<Y^ P(S (2) > x, £_ n > x + ric) + o(F s {x)) 

n 

>x + nc) 

n 

(54) 

+ ^P( J D n , eiR ,^ 2 )>x) + (F s (x)) 
n 

= £i + £ 2 + o(F s (x)). 

Also, 

Si < (1 + o(l))£- £ P(£_ n > x + ric) = (1 + (1))— F s (x)- 

n C 
r<y\ 

On the intersection of the events -D n ,e,_R and {£n. <x — 2R+ (n — l)(c — e)}, 
we have 5"i 2 ^ < x — R. In addition, < R for i = 1, . . . ,re — 1 and, for all 
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i>i 

S^l = S^l + S^-S^<S n+] -S n + S^<R-j{Z-e) + x-R<x. 



(2) 

Thus, on this intersection, Sm < x for all m if x > R. Therefore, 

S (2) >x)= P(D n ^ R , 5 (2) > x, ei 2) > x - 2R + (n - l)(c - e)) 
<P(ei 2) >x-2i?+(n-l)(c-e)). 

Hence, 

£ 2 < £ P(ei 2) > x - 2i? + (n - l)(c - e)) 

n 

= (1 + o(l))—F\x - 2R) = (1 + (1))— F», 

c — e c — e 

as x — ► oo, because of (52). Since e > is arbitrary, the result follows. □ 

Case = . We assume vf = Var(<yW) to be finite for i = 1,2 and 
we use the notation v = \jv\ + v\. 

Theorem 11. Assume a > &W — ft( 2 ) = b and ~P(a^ > x) ~ d^T(x) as 
x — > oo with d^ + d^ > 0, where both F and F s are subexponential. Then, 
as x^ oo, 



P(W^ >x) = 2dW /°° F(x + y(a - &))$ f — ) % 
Jo \v^/yj 

(55) 



+ ^- J F S (r E )+o(F s (x)), 



where is i/ie tai/ o/ the standard normal distribution. In particular, if 
either: 

(i) d^ > or 

(ii) d< 2) = 0, d« > and 

(56) liminf F s (j; 2 )/F s (x) > 0, 

a;— >oo 

then one can replace the equality in (55) by an equivalence and delete the 
term o{F (x)) in this equation. 

Remark 7. Under condition (56), the integral in the right-hand side 
of (55) is of order @(F (x 2 )) = Q(F (x)). Condition (56) is satisfied if the 
tail F (x) is "extremely heavy," for example, F (x) ~ (logx) - ^, K > 0. 
However, if limsup^^ F s (x 2 )/ F s (x) = (the latter holds for Pareto, log- 
normal and Weibull distributions), then the integral is of order o(F (x)). 
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Proof of Theorem 11. We use again (48) and (49), and we are left 
with the problem of finding the asymptotics for 

p(^( 2 )>x,4 l ))= J2 p (^ (2) 

> x,A£l) + o(F s (x)) 

n>N x 

for an appropriate event satisfying the assumptions of Theorem 8. Let 
S m ,n = Y%Li{b - fin+i)- Due to the LLN, maxi< m < n S m , n /n -> in proba- 
bility as n — > oo. Therefore, there exists a nonincreasing sequence rjn^ — > 
such that P(maxi< m < n S mjn /n > r/i 1 ^) < rjn for all n. Take 

KM = { max S m , n /n < rtf*} n {W®^ + ai 2 ^ < x + n(a - b + r^)} 

L l<m<n J 
n,x ' ' n,x 

Easy calculations based on the fact that W^+ m — a - n + Sm,n — (m + l)a 
show that on the event An)c = Kn}x H {c^n > x + n ( a — 6 + ^i 1 "*)}, one has 
Win +m > 0) for all m = 1, 2, . . . ,n. Using the fact that t_2_i > cri 1 ^, one 

(2) 

gets immediately that on this event, = 0. Therefore, on this event, 

rj 2 ^ = oj+i for all j = — n + 1, . . . , and 

W {2) = WP ] = max Sj = V n . 

From the central limit theorem for the reflected random walk, 

weakly, where ip has the following tail distribution: 

P(V> >x) = 2<f>{x) = -= / exp{-y 2 /2} dy. 

V Z7T Ji 

Take any c > 0. If iV x . < n < cx 2 , then 

1 > * > » U(l + o(l))2?f 1 



as x — > oo. For any A > 0, choose c<Cl such that 3>(;^7j) < A. Then 

9 9 

2 P(W (2) > x, A®) < P(K > x)P(ai 1 l > x + n(a - b + V ^)) 

N x N x 

<l + ol rAF s (x). 

a — o 
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If n > ex 2 , then 

V n x \ f x 



P -^=> 



(l + o(l))2$ - 



as x — > oo, since 



V\/n V\lnl \V\/n 



P(-^>y)=(l + (l))2$(y) 



V\ n 



uniformly in y on a compact set. 
Therefore, 



n=cx* 

oo 



J2 P(W^>x,A^) 

E P (^-n-i > x + n(a - 6 + n n ))P(i^ 2 ))P(K > s, tffe 1 )) 

ex 2 

(1 + E F(x + n(a - b + Vn )) ( P f-^= > - o(l) 

CX* 

OO / s 

2(1 + o(l))d« e^( x + n (« - fe )) ¥ ( -4= ) + °(f\ x )) 

2d« F(x + y(a — b))$>( — ^— ) dy + o(F s (x)) 
J ex* \v^Jy) 



poo 

2d (1) / F(x + y{a-b))$ 
Jo 



W X 



v^Jy 



dy 



rf(i) 

+ (F s (x))-A^-0(F s (x)). 
a — 6 

Letting A to 0, we get the result. □ 
Case &« <&( 2 ). 

Theorem 12. Awme < b^ . Then 

(57) P(*W > «) = j^y J"<«) + («£^y ) + 

In particular, if either: 

(i) d^ 2 ) > or 

(ii) d( 2 ) = 0, d?> > and 

(58) liminf F s {2x)/F\x) > 0, 
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then one can replace in (57) the equality by an equivalence and delete the 
term o(F (x)) in the right-hand side. 

Proof. Take the notation from the beginning of Section 6.1 and from 
the proof of Theorem 11. Recall that we consider the case 

6(2) = ft. p u t 

Recall that, from (48) and (49), we have to find the asymptotics for 

n x ,l n x ,2 — l oo 

]r p(w^> x ,Ai i i)= g + £ + x: =pi(x)+p 2 (x)+p 3 (x). 

n>Af;c n=iVa; n=n x , i + l n=n x ^ 

For AT a < n < n x>1 , put K$ x = + < x}. Then, for rfh ] = 

Vn > ~(a - b), on the event a£) x = K$ x n {cr_^ > x + n(a - 6 + %)}, 



Wl 2 2 = and W (2) <K, 



since tj > oj+i. Therefore, 



n x ,l 

Pi(x)< V{a-l>x + n{a-b + ri n ),V n > 



n=N x 



Y, P(^i 1 l>x + n(a-6 + 7 ?n ))P(K> 



n=N x 



^(1 + ,^ , , 

a — b 



since V n /n — > (b — &W) 1 as n — > oo and 

V nx l n xl 

— = ; ■ ► 1 — e < 1 a.s. 

x n x> \ x 

asu oo. 

Consider P2(x). For any sequence rj n — ► and for x sufficiently large, 

"•1,2-1 

P2{x)< J2 P(^l>x + n(a-b + rj n )) 
n=n x i+l 



(l + o(l))dW /-*+n«,2(a-6) _ 

F(x) ax 



a — 6 Ja5+n x ,i(a-6) 

< ( 1+ °y )n ^- n ^r ( x) = (i + o( i))^r ( x), 

a — b n X 2 a — b 



since -F(x) is nonincreasing. 
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Finally, consider P%(x). We will show that, for the appropriate sequences 
{K&% and {rjft}, 

. . tfW - s f fa-b^ e(a-b) 
(59) P 3 (x)~ -—F s l x (— m + 



a-b V \b-bW b-bO) 
Obviously, 

Ps(x)< ?(<T { -l>x + n(a-b + riP)), 

n>n x ,2 

where the right-hand side of the latter inequality is asymptotically equivalent 
to the right-hand side of (59). Now we establish the lower bound. 
From (10), (42), (44) and (45), 

W V=Z-wV-*P-aP 

^ rr c( X ) (!) ( 2 ) 

> ^[-n,o] - max S) - -ay -ay 



. (i) 

> a' — na 



+ max ( ± 

-n<q<0 \ •<— ' 

\m=—n+l / 



0<j< 

Due to the SLLN, as n —> oo, 

1 



max (0, + <7li) - max Sj 

0<?<n— 1 J 



max 5, a) + + <rf ) + ^ ) = r„ 

n n \ 0<j<n-l J n u u 



a.s. 



and 

n -n<9<0 \ ' / 

\m=— n+1 / 

Choose a sequence 5 n J, 0, 77.<5 n — ► oo such that, for all n, 

P(u n >b-5 n ,r n < 6 n , < n(& (1) - a + S n )) > 1 - <5 n 

and denote the latter event by Kn)x (it does not depend on a^l). On this 
event, 

W (2) > oi^ - n(a - b + 2<5 n ) - max (0, + n(b {1) -a + S n )) 
> min (o^ - n{a - b + 25 n ), n(b - - 3<5„)). 
Since n > jz$r> we § et 



n{b — b^ ) — 3n5 n >x-\ — - - — -n — 3n5 n > x 
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for all sufficiently large x. Put rjn^ = 2<5 n . Then, on the event 

= W-l >x + n(a-b + V W)} n tf« , 

we get W > x, and (59) follows. 

Letting e to completes the proof. □ 

6.2. Multiserver queues. The aim of this section is to derive upper and 
lower bounds and sharp asymptotics for the tail of the stationary maximal 
dater of multiserver queues. However, we do not obtain here asymptotics for 
the tail distribution of the stationary waiting time. It is known (see, e.g., 
[17] and [21]) that these asymptotics may, in general, differ significantly 

Since (AA) does not hold, we cannot use the approach of Section 4. We 
show how the ideas of Section 5.2 can be used to derive upper and lower 
bounds which are specific to this queue. 

Recall that we can consider a D/ Gl/m/oo queue with constant interar- 
rival times a. Let Ecr = b and p = G (0, 1). Assume further that P(oi > 
x) = F(x), where both distributions F and F s are subexponential. 

Theorem 13. Under the foregoing assumptions, when x tends to oo, 
P(Z>x) 



a \ma — b a J \b — (m — l)a 

Note that the second term in the right-hand side of (60) disappears when 
b < (in — l)a. 

The proof consists of three steps: 

First, we get a lower bound by using the SLLN. 

Then we get an upper bound by using results from Section 4.3.1. 

Finally, Theorem 7 gives us the tool to derive the exact asymptotics. 

Lower bound. Clearly, 

P(Z>x) > (1 + o(1))p( {J{a- n >x + na}\ 

V n=0 ) 

~ ^2 P(°"-n > x + na) ~ -F s (x). 



a 



Upper bound. Take a sufficiently large L and consider the L-upper- 
bound D/ GJ/l/oo queue with interarrival times La and service times {s n } 
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with mean b = E s\ . Since 

L 

max a, < s i < > <7; , 

\<i<L i 

we get P( ?i > x) ~ LP (a i > x) = LF(x) as x — > oo. Note that, for the mul- 
tiserver queue, 7(0) = b/m. Therefore, we get a natural analogue of Theo- 
rem 5, 

r ^ r r P(R>x) 1 - s 
(blj lim sup g < nm lim sup s = — b [x). 

a:->oo F [xj L^oo x^oo F (x) a — b/m 

Thus, we are in a position to make use of Theorem 7. The rest of the proof 
is quite technical and in the same spirit as that of Theorem 9. Because of 
that, it is omitted. 

APPENDIX 

A.l. Proof of (24). Put b = min (bC* 1 ) ,&(*»)). From Corollary 5, we know 
that, for r = 1, 2, 

{R^>x}c |J (Ag^UB^), 

m r > 1 

where 

A^ = {Y^>x + m r (b^-a)} and P(flW) = o{ F\x)). 
Therefore, 

P(i? (il) >x,R {l2) >x) 
00 00 

< J2 X) P^mi >^ + "ii(fc (il) -a),Y^ >x + m 2 (b (i2) -a)) 

mi=l m2=l 

+ o(F s 0r)) 

(62) < ^ P(F 1 (il) >x + m 1 (b (il) -a))P(y i (i2) >x + m 2 (b (i2) -a)) 

mi^m2 



+ ^ P(min (y i ttlJ ,y i tta >x + m{b- a)) + o( F s (x)) 

m=l 



1 J M 

<^P(y! (il) >i + mi(6 (!l) -a))^P(y i (i2) >x + m 2 (b fe) - a)) 

m 

00 

53 o( F(s + mb)) + o( F s (x)) 



m=l 



Q((F s (x)) 2 ) + o(F s (x)) = o(F s (x)). 
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A. 2. Relaxing the independence assumptions. The aim of this section 
is to give conditions under which assumption (H) of Section 4.1 is satisfied, 
although the r.v.'s Y^> are not independent. 

We assume that there exists a random variable v taking values in an 
arbitrary measurable space (y,By) and such that: 

(i) 

• Given v, the random variables Yf ,j = l,...,r, are conditionally inde- 
pendent. 

• For any j = 1, . . . , r, 

(63) P(y i (i) > x\v) ~ dl j) F(x), 

Py-a.s., where d„ is a nonnegative random variable with a finite mean d^ . 
Then 

(04) ^» ssnp W^M 

x F(x) 

is an a.s. finite random variable, too. 

Assume in addition that, for any 1 < j\ < ji < r, 

32 

(65) E J] d { J ] < oo. 

3=31 

LEMMA 8. Under the foregoing assumptions, for any 1 < j\ < ji < r, 

/ 32 \ / 

.72 

Cj) 

J=Ji 3=51 



pl g y x u; >xj ~p(^ max y^- 1 >z 

(66) 

K ' 32 32 



^p(y 1 (i) >x)~^ ( i (j ' ) F(x). 



Proof. Without loss of generality, we prove the result for j\ = 1, j'2 = r. 
Note that 

P(£i*i 0) >x\u) ^ {j) ^ P^axj-yfi >x\u) 
F(x) ~T " ^~ F(x) 

P^-a.s. and, for all x, 

P(ma Xi ytt> > xH < P(Ei ^ > 4^) < A Zi) mn TM 
- F{x) ~ F(x) ~\ [ u x P F(x) ' 
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where the latter supremum is finite. Then the dominated convergence theo- 
rem implies that 

and 

P(maxjY^>x) 

W(x) ' □ 

Consider the following example, which covers the generalized Jackson 
network case. Assume that there are given: 

• Some random vector v= (i/W , . . . , iM> ) with nonnegative integer- valued 
components, such that Eexp(cf^-') < oo for some c > and for all j = 
l,...,r; 

• r sequences {<Jn } of i.i.d. subexponential random variables that are mu- 

ff) 

tually independent and independent of v, and such that ~P(a\ > x) ~ 
fo>F(x). We do not make the assumption that the r.v.'s v^ l \ . . . , are 
independent. 

Put = J2i=i a i- The above conditions imply that 

Eexp ^-yV^ < Eexp ^cmax^'j < ^Eexp(czz^) < oo 

and that for all j = 1, . . . , r, 

(j) p (°i ] > <) 
u u ' = sup = < oo. 

t F(t) 

Due to subexponentiality, for j = 1, . . . , r, 

It is known (see, e.g., [13], page 41) that, for any e > 0, one can choose 
K® = KW(e) such that 

P(y i C?) > x\u) < K {j \l + eY U) P(a{ j) > x). 

The right-hand side of the latter inequality is not bigger than K^u^\l + 
e) uU) F{x). 

Take e > such that log(l + e) < d . Then the conditions of Lemma 8 are 
satisfied with d { J } = v®lV\ d^ = 0EvW and dV } = K^u^\l + ef U) . 
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A. 3. Deterministic interarrival times. We extend to the monotone sep- 
arable framework the approach used in [4] for single server queues to show 
that there may be no loss of generality in assuming that a network has deter- 
ministic interarrival times when one wants to evaluate the tail asymptotics 
of its stationary maximal dater. 

The framework is that of Section 2. Fix {Cn}, {fi} and consider a family 
of networks with different "input sequences" {T n } such that En > 7(0). 
Without loss of generality assume To = 0. 

In particular, a network with constant interarrival times (say a) belongs 

to this family. For such a network, we use the notation and Z^_ n j . 
For any {T n } and for any a < Eti, set 

ip({T n }, a) = sup(na + T_ n ) = sup V (a - n) . 

n>0 n>0Vi=_ n / 

Theorem 14. Assume that there exist a continuous and strictly posi- 
tive function h : (7(0), oo) — > (0, 00) and a subexponential distribution G such 
that, for any a > 7(0), 

(67) P(Z (a) > x) ~ h(a)G(x) as x — > 00. 

Then, for any network with random interarrival times {r n }, such that Eri = 
a > 7(0), the following is valid: If {r n } and {Cn} ore independent and if, for 
any a < a, 

(68) P(V>({r n },a) > x) = o(G{x)) asi^oo, 
then 

(69) P(Z>x)~/i(a)G(x) asx^oo. 



Remark 8. In particular, condition (68) is satisfied if the r n 's are i.i.d. 
Indeed, then tp({T n },a) has either abounded [if P(5 > t\) = 0] or exponential 
tail, which is lighter than any long tail. 

Proof of Theorem 14. Take any e £ (0, a — 7(0)). Due to the mono- 
tonicity, 

Z[~n,0] < Zfc* + max^ ^(a-e- n ) 

Therefore, 

Z < Z^ + ^({r n },a -e) = Z^ + ^ 
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where Z^ a e ' and ip are independent. Therefore, 

P(Z >x)< P(Z^-^ + 1\> > x) ~ P(Z( a " £ ) > x) ~ /i(a - e)G(x). 

Thus 

hm sup — =— — < h(a - e) 

for any e G (0, a — 7(0)). Letting e go to 0, we get the upper bound /i(a). 

For the lower bound, we use the monotonicity, the SLLN for the r's, the 
LT and the independence assumptions. For any e > 0, one can choose a 
sufficiently large C = C{e) such that 

P(T_ n > -n{a + e) - C V?i > 0) > 1 - e. 

Denote the latter event by D e . Then 

P(Z >x)> P(Z > x, D £ ) > P(Z( a+e) - C > x, D £ ) 

> p(z (a+£) - C> x)(l - e) ~ h(a + e)(l - e)G(x + C) 

~h(a + e)(l-e)G(z). 

Thus, for any e £ (0, 1), 

PC z > rl 
l im in f r \*>*) >fe(a + e)(l-g). 

Letting e go to 0, we get the lower bound with coincides with the upper one. 
□ 
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